Anisotropic flow in 4. 2 A GeV/c C+Ta coflisions 



Lj. Simic ^ and J. Milosevic 
^ Institute of Physics, P.O. Box 68, 11081 Belgrade, Yugoslavia 
^ Physikalisches Institut, Philosophenweg 12, 69115 Heidelberg, Germany 
(Dated: February 8, 2008) 

Anisotropic flow of protons and negative pions in 4.2A GeV/c C+Ta collisions is studied using 
the Fourier analysis of azimuthal distributions. The protons exhibit pronounced directed flow. 
Directed flow of pions is positive in the entire rapidity interval and indicates that the pions are 
preferentially emitted in the reaction plane from the target to the projectile. The elliptic flow of 
protons and negative pions is close to zero. Comparison with the quark-gluon-string model (QGSM) 
and relativistic transport model (ART 1.0) show that they both yield a flow signature similar to the 
experimental data. 



The anisotropic transverse flow of particles has been 
actively studied in nuclear collisions over a wide range of 
energies. At lower energies ||, ||, ||, || , the flow is usu- 
ally studied in terms of the mean in-plane component of 
transverse momentum at a given rapidity, {p^{y)) Q, and 
additionally quantified in terms of derivative at midrapid- 
ity Fy = d{p^)/dy. At high energies, the Fourier expan- 
sion of the azimuthal distribution of particles constructed 
with respect to the reaction plane is used ||, 0, |§1 . In this 
expansion the first harmonic, vi, quantifies the directed 
flow while the second harmonic, V2, quantifies the ellip- 
tic flow. Using the Fourier expansion, the anisotropic 
transverse flow was analyzed for heavy symmetric sys- 
tems at the AGS |, |, 0, SPS 0, |l| and RHIC |l| 
energies. It was found that the flow observables are im- 
portant tools for investigating properties of high density 
region created during the initial collisions. In particular, 
the elliptic flow measurements may provide an important 
constraint on the Equation of State (EOS) of high density 
nuclear matter |l^ . 

In this paper the anisotropic transverse flow of pro- 
tons and negative pions in 4.2A GeV/c C-|-Ta collisions 
is studied using the Fourier analysis of azimuthal dis- 
tributions. The analysis is performed using 1000 C-|-Ta 
semicentral and central collisions obtained with the 2-m 
propane bubble chamber, exposed at JINR, Dubna syn- 
chrophasotron. The semicentral and central collisions are 
selected by rejecting « 50% events with the smallest mul- 
tiplicity of participant protons. Additionally, the same 
type of analysis is performed using 100000 events gen- 
erated by the Quark-Gluon-String Model (QGSM) |l|] 
, and the same number of events generated by the rela- 
tivistic transport model (ART 1.0) [|l^. For these events 
the same centrality criterion is applied as in experiment, 
leading to the average impact parameter (b) « 4.54 (4.05) 
fm according to QGSM (ART 1.0). 

In order to study the inelastic interactions with tanta- 
lum nucleus, (^^^ Ta), three tantalum foils (1 mm thick 
and 93 mm apart) were placed inside the chamber work- 
ing in the 1.5-T magnetic field. The characteristics of 
the chamber allow precise determination of the multi- 
plicity and momentum of all charged particles, as well as 



identification of all negative and positive particles with 
momenta less than 0.5 GeV/c. All recorded negative 
particles, except the identified electrons, are taken to be 
7T~ . Among them remains admixture of unidentified fast 
electrons (< 5%) and negative strange particles (< 1%). 
All positive particles with momenta less than 0.5 GeV/c 
are classified either as protons or tt^ mesons according 
to their ionization density and range. Positive parti- 
cles above 0.5 GeV/c are taken to be protons, and be- 
cause of this, the admixture of 7r+ of approximately 7% is 
subtracted statistically using the number of tt~ mesons 
with p > 0.5 GeV/c as follows: Hp — n+ — nT^+{p < 
0.5 GeV/c) - 0.82 • n^-ip > 0.5 GeV/c), where n+ de- 
notes the number of single positively charged particles, 
and 0.82 takes into account the proton deficit in tantalum 
nuclei and consequently also tt"*" deficit. From the ratio 
for each momentum interval we determine the weight of 
protons which we further use when calculating distribu- 
tions of other kinematical variables. From the result- 
ing number of protons, the projectile spectators (protons 
with momenta p > 3 GeV/c and emission angle 6* < 4°) 
and target spectators (protons with momenta p < 0.3 
GeV/c) are further subtracted. The resulting number of 
participant protons still contains some 17% of deuterons 
(with p > 0.48 GeV/c) and 11% of tritons (with p > 
0.65 GeV/c). The experimental data are also corrected 
to the loss of particles emitted at small angles relative to 
the optical axes of chamber and to the loss of particles 
absorbed by the tantalum plates. The aim of this cor- 
rection is to obtain isotropic distribution in azimuthal 
angle and smooth distribution in emission angle (both 
measured with respect to the direction of the incoming 
projectile). 

The azimuthal distribution of particles may be repre- 
sented with the first three terms of the corresponding 
Fourier expansion 

^ ~ ^ [1 + 2i^i cos((^) + 2v2 cos{2<P)] , (1) 

where the two coefficients, vi and V2, quantify the di- 
rected and elliptic flow via vi — {cos{(/))) and V2 — 
(cos{2(j))). In Eq. (1), </> = tpiab - ^piane IS the parti- 
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cle azimuthal angle determined with respect to the re- 
action plane, with (piab denoting the azimuthal angle of 
particle in the laboratory frame and ^piane denoting the 
azimuthal angle of the (true) reaction plane. Since both 
the projectile momentum and the impact parameter vec- 
tors are available in the QGSM/ART simulations, they 
are used to determine the corresponding reaction plane. 
In the experiment the reaction plane is determined, for 
each event, using the projectile momentum vector and 
the vector Q determined from 

Q = X! P'^'^y > ycm + S) PTjiy < Vcm - S), (2) 

where represents the transverse momentum of the 
proton emitted in the forward (y > t/cm + <5), or back- 
ward (y < Ucm — <5), hemisphere. Here, j/cm denotes the 
center-of-mass rapidity of participant protons while the 
quantity 6 (=0.2) removes the protons emitted around 
the 2/cm which are not contributing to the determination 
of the reaction plane. The reaction plane angle for a pro- 
ton is determined using this expression only if this proton 
is not included in the above sum (i.e. if its rapidity lies 
in the interval from ?/cm — 5 io ycm + 5). Otherwise, in 
order to avoid autocorrelation (which is an effect of the 
finite multiplicity), the Q vector is constructed by the 
analogous expression in which the contribution of this 
proton is simply omitted We found that the reac- 
tion plane angle distribution is essentially flat, thus con- 
firming the absence of significant distortions which could 
influence the magnitude of the extracted flow parame- 
ters. The accuracy with which the reaction plane angle 
is determined, i.e. the reaction plane resolution, is eval- 
uated by the subevent method In this method, 
each event is divided randomly into two subevents, and 
then the corresponding two reaction planes are deter- 
mined. Subsequently, the absolute value of the relative 
azimuthal angle, $12, between these two estimated reac- 
tion planes is obtained. The relative azimuthal angle dis- 
tribution is the basis for the correction of the Fourier co- 
efficient, v'l^ obtained with the estimated reaction plane. 
The relationship between the v'l, and the Fourier coeffi- 
cient vi obtained relative to the true reaction plane, is 
v'l = vi (cos(A$)), where (cos(A$))~^ is the correction 
factor determined from $12 distribution following the 
prescription given in |||, [l^. We find (cos(A$))=0.59. 
The correctness of this procedure is checked using the 
QGSM. Using this model, the coefficient vi vs. rapidity, 
for protons and negative pions, is calculated with respect 
to the true reaction plane and also with respect to the 
estimated reaction plane. The result of comparison is 
presented in Fig.l. 

The QGS calculations show that the vi values obtained 
with respect to the estimated reaction plane, after apply- 
ing correction procedure, are somewhat underestimated 
around projectile rapidity for protons and negative pi- 
ons. The QGS calculations also show that the correction 
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FIG. 1: Rapidity dependence of v\ and V2 for protons and 7r~ 
for 4.2A GeV/c C-I-Ta collisions: top- filled circles represent 
the experimental results for vi while the solid (dashed) line 
represents the QGSM calculation for vi with respect to the 
true (estimated) reaction plane; bottom- filled circles repre- 
sent uncorrected experimental V2 values (see text), while the 
solid (dashed) fine represents the QGSM calculation for 112 
(U2) with respect to the true (estimated) reaction plane. 



procedure for V2 as outlined above is not applicable be- 
cause of the smallness of the elliptic flow. Therefore, in 
the following analysis, this coefficient is not corrected to 
the reaction plane resolution. 

Fig. 1 (top) displays the experimentally determined 
vi coefficient vs. y (with y calculated in the lab frame), 
for protons and negative pions. In the case of protons 
the dependence of vi on rapidity is characterized by a 
curve with a positive slope and with the zero-crossing at 
y « 0.5, that corresponds to average rapidity of protons. 
The curve indicates a positive directed flow with magni- 
tude vi ~ 0.2, at rapidities close to the projectile rapidity 
{yp = 2.2, at p=4.2A GeV/c). The QGSM reproduces 
satisfactorily the shape of viijj) curve and within error 
bars reproduces the magnitude of the flow. The exper- 
imental results are also compared with the relativistic 
transport model ART 1.0. These are shown in Fig. 2 
where the calculations are performed both for stiff and 
soft EOS. ART model yields a directed flow which flows 
trend similar to the experimental data, but underesti- 
mates the flow intensity in the projectile/target rapidity 
region. 

Using the extracted values of vi and their relation to 
the mean transverse momentum projected onto the re- 
action plane, vi = (px)/{pt), we can evaluate (px) as 
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FIG. 2: Experimental results for vi, V2 as a function of 
rapidity, compared with ART 1.0 model calculations. 



a function of rapidity and determine the slope, F — 
dijpx) / d{y / Hh) , with respect to rapidity normalized to 
projectile rapidity in the lab frame. In the present anal- 
ysis we find for the slope at mid-rapidity F — 215 ± 32 
MeV/c. Comparison with the other results obtained 
at the same energy, for various C-nucleus combinations 
shows increasing of the slope with the target mass: F ~ 
144 MeV/c for CC |l|l 134 MeV/c for CNe [|o) and 
198 MeV/c for CCu ||2(|. After the normalization to 
the mass number of the colliding system we obtain the 



so-called scaled flow F; 
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27 ± 4 



MeV/c, that allows a comparison of the energy depen- 
dence of flow values for different projectile/target mass 
combinations. This value is in agreement with the ob- 
served trend |^ that after reaching the maximum at a 
beam energy around 0.7-2A GeV, the directed flow slowly 
decreases with increasing beam energy. 

For negative pions the experimental values of vi are 
positive in the entire rapidity interval. The vi is largest 
in the target rapidity region (u™"^ k, 0.10) and monoton- 
ically decreases with increasing rapidity towards the pro- 
jectile rapidity. Such vi dependence on y reflects the fact 
that the pions are preferentially emitted in the reaction 
plane from the target to the projectile. This behavior 
is attributed to a shadowing effect of the heavy target. 
Both QGSM and ART model cannot strictly account for 
the vi{y) dependence for negative pions. 

Figure 1 (bottom) displays the experimentally deter- 
mined V2 coefficient versus y for protons and negative 
pions. The uncorrected values of v'2 show that in the en- 
tire rapidity interval the elliptic flow is small 1 < 0.02 
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FIG. 3: Rapidity dependence of v\ for protons and tt" orig- 
inating from: primary non- resonant interactions (stars), de- 
cay of resonances {full circles), and secondary non-resonant 
interactions {open circles), for 4.2A GeV/c C-f-Ta collisions 
generated with the QGSM. 



if not zero, and this is consistent with QGSM (Fig. 1) 
and ART (Fig. 2) predictions. The values of v'2 addition- 
ally confirm the result ||lj] that at beam energy of « 4 
GeV the elliptic flow exhibits a transition from negative 
(out-of-plane) to positive (in-plane) eUiptic flow. 

Since the QGSM predictions are in fair agreement 
with the various experimental results at 4.2A GeV/c, 
we use this model to clarify the question of which of 
the processes are responsible for the flow effect. In 
this model, for C-|-Ta collisions « 43% of protons and 
« 83% of vr^ originate from decay of the lowest-lying 
resonances (A's, g, lo, 77 and r\ ). The rest originates from 
the 'non-resonant' primary and secondary interactions 
of the type: NN NNtt, AN AiV, ttN ttN, 
t:NN — > NN. The protons and pions from primary 
interaction escape the collision zone without further 
rescattering and comprise w 1% of the total. According 
to QGSM, we separately evaluate the flow of protons 
and pions originating from the following sources: 
(z) decay of resonances; 
(ii ) primary non-resonant interactions; 
{Hi ) and secondary non-resonant interactions. 

Figure. 3 shows Vi vs. rapidity for protons and nega- 
tive pions originating from the decay of resonances, and 
from primary and secondary non-resonant interactions. 
The protons originating both from the decay of reso- 
nances and from the secondary non-resonant interactions 
show a directed flow of similar intensity. The same ap- 
plies to the flow of pions. For both protons and pions 
from the primary interactions the directed flow has a 
maximum around y ~ 0.6, and decreases towards the 
projectile and target rapidities. These protons and pions 
are produced at the early stage of the collision, and both 
are shadowed by the cold spectators. Later, after the 
spectator matter leaves the collision zone, rescattering of 
protons near the beam (target) rapidity region is small. 
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while the pions are still affected with the shadowing effect 
of the participant nucleons trough both pion rescattering 
and reabsorptions. This could be the underlying mecha- 
nism that leads to the different intensity and dependence 
on rapidity for the directed flow of pions and protons. 

In summary, the directed and elliptic flow of protons 
and negative pions in 4.2A GeV/c C+Ta collisions 
was examined using the Fourier analysis of azimuthal 
distributions of experimental events and also by using 
the events generated by the QGSM and ART 1.0. The 
protons exhibit strong directed flow with magnitude 
vi ~ 0.2 at rapidities close to the projectile rapidity. 
The directed flow of pions is positive in the entire 
rapidity and slightly peaked at target rapidity, where 
vi w 0.1. This behavior indicates that the pions are 
preferentially emitted in the reaction plane from the 
target to the projectile. For both sets of particles in the 
entire rapidity interval the elliptic flow is close to zero 
{\v2\ 0.02), this being consistent with the result that 



at b eam energy of rj 4 GcV the elliptic flow shows a 
transition from negative to positive. A comparison with 
the quark-gluon-string model (QGSM) and relativistic 
transport model (ART 1.0) shows that they both yields 
the flow signature similar to the experimental data. 
Additionally, the QGSM shows that two factors that 
dominantly determine the proton and pion flow at this 
energy, are the decay of resonances and the rescattering 
of secondaries. The shadowing by the cold spectator 
matter affects only the flow of the particles produced at 
the early stage of the collision. 

The authors are grateful to members of the JINR 
Dubna group that participated in data processing. 
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